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§ ■ ABSTRACT 

(N : 

I Many hydrodynamic processes can be studied in a way that is scalable over a vastly 

I^H I relevant physical parameter space. We systematically examine this scalability, which 

' has so far only briefly discussed in astrophysical literature. We show how the scalability 

,— i' is limited by various constraints imposed by physical processes and initial conditions. 

' Using supernova remnants in different environments and evolutionary phases as appli- 

O: cation examples, we demonstrate the use of the scaling as a powerful tool to explore 

. the interdependence among relevant parameters, based on a minimum set of simula- 
tions. In particular, we devise a scaling scheme that can be used to adaptively generate 

2 I numerous seed remnants and plant them into 3D hydrodynamic simulations of the 

' supernova-dominated interstellar medium. 

^ ■ Subject headings: methods: miscellaneous — galaxies: ISM — ISM: structure — super- 

^ ' nova remnants 

m : 
o . 

^ : 

, 1. Introduction 
<N ■ 

■ Similar natural phenomena, which may occur on vastly different space and time scales, can 

■ often be treated in the same way. One well-known example is the similarity between a supernova 
1^ . (SN) in the interstellar medium (ISM) and a nuclear explosion in the earth atmosphere despite of 

their vastly different energies (~ 10^^ ergs vs. ~ 10^^ ergs). Their blastwave structure and evolutio n 



X 



may be mathematically approximated by the same self-similar Sedov- Taylor solution (jSedov Ill959l ). 
I with an appropriate scaling according to the energy and ambient medium density. Such self- 

similarity, though often limited in its applicability (e.g., the evolution needs to be adiabatic; the 
mass of the ejecta is negligible; etc.), has been widely used in astrophysical studies. 

The scalability of a hydrodynamic process, as will be demonstrated in the present paper, has 
a much broader application. Here we explore how the solution (or simulation) for one physical 
setup can be scaled to another, when the underlying governing equations are the same. The self- 
similarity is then only a special case of the scalability. Therefore, the scalability analysis provides 
a systematic way to examine the physical parameter space, based on a limited number of solutions. 
As a specific example, we apply our scalability analysis to the study of the SN remnant (SNR) 
evolution in various environments and at different evolutionary stages. 
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The scalability has the same idea as the homology relations, which are used in studying the inte- 

rior s tructure of stars in complete equilibrium (both hydrostatic and thermal; e.g.. lKippenhahn &: Weigert 
lOQ^l . A unique s c aling relation to study the SNR evolution was probably first introduced bv lSgro 
(|l972l ). IChevalierl (jl974l ) discussed the same scaling relation to analyze the evolution of SNRs of 
different setups with a limited number of simulations. By recognizing that one simulation of a 
particular SNR can be used to describe a family of SNRs if they all have the same Egnn^ (where reo 
is the number density of ambient medium and Egn is the SN energy, see §2.2 for further discussion), 
Shelton et al. I (jl999l ) pointed out the usefulness of the scaling in interpreting observations with a 
few simulations. These discussions, though limited in their scope, have demonstrated the potential 
of using the scalability in the study of SNRs. 

In the present paper, we attempt to give a systematic examination of the scalability of SNR 
solutions and simulations and provide specific application examples. The initial motivation of this 
work is to find an effective method to generate ID SNR seeds that can be embedded into 3D 
hydrodynamic simulations of the SN-dominated ISM, particularly in galactic bulges where the ISM 
is dominated by diffuse hot gas. The SNR evolution in such ambient medium in general cannot be 
described by the self-similar Sedov- Taylor solution, which assumes a cool ambient medium (hence 
with no energy co ntent). In fact the e volution depends on both the density and temperature of the 
ambient medium ( Tang &: Wane]|2005l ) . Each 3D simulation needs, for example, more than 10^ SNR 
seeds for a bulge of an even moderate stellar mass ~ 10"^" as in our Galaxy or M31, over a few 
times their dynamic time scales (~ 10^ years). The seed embedding, worked with an adaptive mesh 
refinement scheme, can effectively extend such a 3D simulation to include the subgrid evolution of 
SNRs. Here the subgrid evolution means that the structure of SNR seeds results from the evolution 
on scales much smaller than the highest spatial resolution available in the 3D simulations. The size 
of an embedded SNR seed cannot be too big (in order to use the ID simulation properly) or too 
small (to be within the limited dynamic range of a 3D simulation). Therefore, we should adaptively 
select suitable SNR seeds according to the local density and temperature (values and gradients) 
of the environments. For each selected SNR seed, we need the ID radial density, temperature, 
and velocity profiles, with proper normalizations to gu arantee the mass, momentum, and energy 
conservations of the embedding into the 3D simulation (|Tang et al. II2009I ). In principle, we could 
draw the seeds (with some interpolations) from a library of the profiles in a grid of the three 
parameters: SNR radius as well as the density and temperature of the ambient gas (the explosion 
energy and ejecta mass are assumed to be the same for all SNe; otherwise a larger parameter space 
is required for such a library). Clearly, this approach of generating and using such a large library 
is not elegant, if practical. Instead, we find that we can use the scalability to generate the seeds 
based on a very limited number of ID SNR simulations. We describe this simple approach as an 
application example. 

In §2 we show how to deduce the scaling relation starting from the basic gas dynamics equations 
and what the constraints of the scaling are. In §3 we apply the scaling to specific cases of the SNR 
evolution. In particular, we demonstrate how we use the scaling to generate SNR seeds for the 
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3D simulations with a few ID simulations and how to correctly interpret the simulated relations. 
Finally in §4, we discuss the potential use of the scaling in a broader context. 



2. Scaling Scheme 
2.1. Basic Idea 

For a system passively evolving without source terms, the dynamics can be described by the 
following equations: 

|^ + V.(pv)=0, (1) 

^ + V.(pvv) + VP = 0, (2) 

^ + V . [{pe + P)v] = 0, (3) 

where p, v, and P denote density, velocity vector, and pressure, while the total specific energy e 
can be expressed as e= (^^i)p + for ideal gas. This set of equations is in a closed form; i.e., 
we can solve five independent unknown scalar variables {p, v, and P; t and the spatial position 
are explicitly known variables), from the five equations (because Eq. 2 can be decomposed in three 
scalar equations). If the thermal state of the gas is of interest, the equation of state is needed: 

P = R,.PT (4) 

where is the equivalent ideal gas constant^ and T is the gas temperature. The scaling scheme 
can be found by converting each variabl^ Q as QX'^ (where A can be any positive value) and 

solving them for iq (to recover the equations before the conversion). It can be found that the same 
solution holds if 

ip-it = ip + iv-iL, (5) 

ip + iv - it = ip + 2«t, -ii = ip- ii, (6) 

ip + ie - it = ip + ie + iv - ii = ip + iv - ii, (7) 

ip = ip + iT, (8) 



^Note that i?^ — k/firUp, where k is the Bohzmann constant, rup is the proton mass, and fi is the average atomic 
weight. The value of fi depends on the gas ionization state and might change with temperature. For hot gas n only 
weakly depends on temperature so Rp, can be approximated as a constant. And as long as the scaling is within a 
limited temperature range, the small variation in fi can be neglected. 

^These variables include fundamental quantities such as length (L), time (t), and mass (M) and other physical 
quantities such as density (p), pressure (p), velocity (v), specific energy (e), total energy {E), etc. 
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which can be simphfied as 

iv = iL- k, (9) 

ip = ip-2iL + 2it, (10) 

ie = 2iL-2it, (11) 

iT = 2iL-2it. (12) 

Therefore, only three indices are independent in Eqs. ([9]) — p^ . If we consider ii, it, and im 
to be independent, other indices can then be expressed as 

iv = iL- it, (13) 

ie = 2iL - 2it, (14) 

V = Uf - 3iL, (15) 

ip = iiM - ii - 2it, (16) 

iE = iM + 2iL-2it, (17) 

ir = 2iL - 2it. (18) 

Eqs. (jl3p - (jl7p show that the scaling relation of the physical quantities can be directly inferred from 
their dimensions based on the basic units of mass, length, and time. Indeed there is no constraint 
on the choice of if , it, or ijyj from the governing equations (l)-(3) for this simple case (see also 



Rvutov et al.lll999l where such a property of the equations is called Euler similarity). But ix is 
restricted by and it through the equation of state. This constraint is given because is fixed 
(to a number with non- vanishing dimension hard-wired in a specific simulation), which reduces 
one degree of freedom for the scaling of the pressure, density, and temperature. In other words, 
although we have four basic units for the ideal gas hydrodynamics (i.e., mass, length, time, and 
temperature), we are only able to freely change three of them when scaling from one case to another, 
i.e., za/ and two other indices from the pool of iL, it, and ix- 

A special class of the scalability is the self-similar solution. In this case, clearly only one 
solution is needed. However, such a solution, if exists, may not be easily expressed in an analytic 
form and may be applicable only asymptotically (e.g., when the effect of the initial condition 
becomes negligible). In general, one may resort to a simulation to reach the solution. Thus it can 
be studied as part of the scalability problem considered here. 



2.2. Additional Constraints 

If Eqs. (l)-(3) have source terms, more constraints may then be placed on the scaling relation. 
For example, the inclusion of the thermal conduction term, g = V • [k{T)'VT], at the r.h.s in 
equation (3) requires 

- 3if -I- - 3.5zT = (19) 
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for non-saturated thermal conduction in which fi:(T) = koT^^'^, where fcp ^ 9x 10 ^ergcm 

is the Spitzer conduction coefficient (jSpitzer Ill962l ). For saturated thermal conduction {q oc pc^) 



no constraint like Eq. (jl9p is required because it does not require an extra coefficient with a 
non-vanishing dimension. The constraint from a radiative cooling term — n^neA(T) in the same 
equation depends on the form of emissivity A(r). For optically thin primordial gas of temperature 
larger than 5 x 10^ K, for example, the emissivity can be approximated as A(T) = AqT^/^ where 
Aq ~ 10~^^ erg cm^ s~^ K~^/^ and the constraint becomes 

^M + 3i^-5^L + 0.5^^ = 0. (20) 

In general the cooling rate A(T) does not have such a simple power law form, so the scaling relation 
is 

ir = 0; iL = it; = 2it (21) 



when combined with Eq. (fT8]) , which gives the unique scaling relation adopte d bvlSgrd (119721). Th is 



scaling relation also makes EstiTIq an invariant (i.e., ie + '^ip = 0) as used in lShelton et al. I ( 19991 ) 



Additional physical constraints other than those from the governing equations may need to be 
placed on the scalability. For example, the scaling requires iE=^ and/or ip=0 between solutions 
with an identical explosion energy and/or ambient density. Note that we have three degrees of 
freedom for all the power indices, if the number of the constraints is less than three (i.e., at least 
one index is free to change), the solution is then scalable; otherwise the solution pertains only to a 
particular problem. 

Implicit constraints on the scaling relation may be imposed by initial conditions as well. Specif- 
ically, when we scale one solution [pa{ra-,ta), Ta{ra,ta), ■■■] to another [pb{rb,tb), Ti,{ri,,tb), ■■■], the 
corresponding initial condition needs to be scaled in the same way. For example, a particular scaling 
relation can be determined by specifying im, ii, and it, which in turn determines i^, ip and other 
indices via Eqs. (fT3]) - (fT8| ). This scaling relation demands that the corresponding initial conditions 
should be related by 

Pb{rb,tbo) = Pa{ra,tao)^''', (22) 

Vbin^tbo) = Va{ra,tao)y" , (23) 



and other quantities for 



r6 = r,A^^ tbo = taoX''. (24) 



It is such demands on the initial condition that often make one problem be unique from others 
(limiting the scalability of their solutions), even if all have the same governing equations and 
characteristic quantities such as total energy and mass (see §3.5 for further discussion). In the 
following we assume that the scalable solutions do have the required initial conditions unless being 
explicitly expressed otherwise. 
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3. Application Examples 

We use the evolution of SNRs as a simple example to demonstrate how the above described 
scalability can be used. The scalability of an SNR solution or simulation depends on its evolutionary 
stage and on the properties of the ambient medium. 



3.1. Sedov- Taylor Solution 

If the ejecta mass can be neglected and the ambient gas temperature can be approximated to be 
zero, then the evolution of the SNR can be described by the Sedov- Taylor solution, which depends 
only on the explosion energy Egn and the amb ient gas den sity po- The solution can be obtained 



either numerically (Taylor 1950) or analytically (jSedov Ill959l ). In particular the self-similar solution 
of the shock front, 

rsH{t) = i J , (25) 

is widely used, where ^ ~ 1.15 for ideal gas with the specific heat ratio 7=5/3. From the scaling 
point view, following Eq. (jl5p and (jl7p . we have fyiL = ie + "^^1 — ip (i-e., r oc 

^1/5^2/5^-1/5)^ which 

just the same relation shown in Eq. ()25p . Furthermore, for the same remnant, we have i_B=0 and 
ip=0, hence ZL=0.4if, i„=— 0.6z(, ip=—1.2it, and other indices following Eqs. (fT3]) - (fT8]) . It shows 
that for a self-similar solution all the non-zero indices are proportional to it . This allows the scaling 
from one solution at any particular time to another. 

Of course, the Sedov- Taylor solution applies only when the SN ejecta and ambient temperature 
can be neglected. Otherwise, this self-similar solution cannot be applied. But the scaling may still 
be useful. 



3.2. SNRs in Hot Gas 



If the ambient t emperature is not n egligible, a generalized formula for the SNR shock front 
can be expressed as (jTang k. Wangll2005l ) 



rsh{t) 



E 

Po 



2 \ 1/5 



F 



3 2 7 t 
"5' 5' 5'~t^ 



(26) 
(27) 



where Cg is the sound speed of the ambient medium, F is the generalized hyper-geometric function 
and is equal to 1.16 when t = tc which is defined as 



2 A Esn 
5 J pocf 



1/3 



(28) 
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This modification accounts for the energy content of the swept-up ambient medium. Note that 
Eq. (j27p is very similar to the Sedov- Taylor solution Eq. (|25p except for the modification term F. 
When t > tc, the shock front evolution significantly deviates from the Sedov- Taylor solution. If 
the temperature of the ambient medium is zero, then tc ^ oo, F = 1, and Eq. (j27p is the same 
as Eq. (|25p . In general, the solution in this case is no longer self-similar. The reason is that the 
evolution requires i_E=0, ip=0, and ip=0. Thus all the indices are fixed to be zero. The internal 
profiles change with time and cannot be scaled from one time to another. 

But the solution is still scalable between remnants evolving in different environment. As 
illustrated in the following, one simulation in a particular environment is sufficient to infer specific 
SNR solutions in other environments with different ambient density and/or temperature. 



3.3. SNR Seed Generation for 3D Simulations 

The scalability has a particularly important application in planting SNR seeds in 3D simula- 
tions of the ISM. Suppose that each SNR in such a 3D simulation has the same explosion energy 
(i.e., iE= For convenience, we can choose the remaining two free indices to be the power in- 
dices of the density and temperature, which can be directly measured in the simulation. The scaling 
relation can then be simplified as 

iE = 0, (29) 

IM = -IT, (30) 

iL = -it/^ - ip/3, (31) 
it = -5zt/6 - ip/3, (32) 
= O.bix- (33) 

Therefore, we can build a library of SNR templates. Each consists of the radial profiles of density, 
temperature, and velocity when the shock front of the SNR has a certain radius or age. These 
templates can be obtained from a ID simulation of an SNR evolving in a uniform ambient medium 
of density pa and temperature T^. Using the library and the above scaling relation, we can generate 
SNR seeds at any time and at any position of the 3D simulation. The time and position of each 
SN can be realized randomly according to the Poisson statistics and the stellar distribution of a 
galactic bulge, for example. The procedure to embed an SNR seed into the 3D simulation is as 
follows: 

1) At the time step just after the SN and around its position, determine a spherical region of 
radius r^ax , within which T and p are sufficiently uniform so that a ID SNR seed is a reasonable 



' In principle, the explosion energy can vary as well and the resultant scaling relation can be obtained in a similar 
way, and we do not need to expand the parameter space of the SNR library. Without losing generality, however, we 
have assumed the canonical value Esn ~ W''^ erg for Type la SNe. 
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approximation (in practice, a fraction of the radius r^, = rjrmax niay be used, where rj < 1); 

2) Calculate the average density p and (mass-weighted) temperature T in that region, and then 
determine the it and ip-. ip = logx{p/pa) and ir = logx{T/Ta), which then determines il (Eq. [3T]) : 

3) Search in the library for an SNR template which has the shock front radius of ~ r^A"*^ 
and the corresponding SNR age t^; 

4) Wait to a future elapsing time tb of the simulation, when tb — iaA*' is just satisfied, and read 
from the library the template (which may be interpolated to account for the difference between 
tb and taX^\ though not necessary if the time step is sufficiently small, compared to the tb); 

5) Scale each profile Q of the template according to QA*<3; 

6) Plant the scaled profiles int o the 3D simulation by replacing the values within of the SN 



position (see lTang et al. II2009I for more details). 



As the result, we have an adaptively configured SNR seed in the 3D simulation. This seed has a 
dynamically self-consistent structure expected for the SNR evolving in the local ambient medium, 
which is particularly important for accurately tracing the SNR structures and SN ejecta (see §3.4). 
We can therefore incorporate the sub-grid evolution of the SNR into the large-scale 3D simulation, 
which optimizes t he use of the computational time and enlarges the covered dynamical range 



(jTang et al. II2009I ). It is not clear how such realism of the SNR seed and the adaptiveness of its 
planting can be realized in other simple way (e.g., assuming a uniform thermal energy deposition 
or other arbitrary profiles). 



3.4. SN Ejecta and Scalable Initial Condition 

When an SNR is young, the mass of the SN ejecta can be considerable. Assuming that the 
mass is the same for the SNRs in the consideration, we have iM=0 as well as i_B=0 and ip=0, 
as in the previous case. The SNR evolution is not self-similar and asymptotically approaches the 
Sedov- Taylor solution only when the swept-up mass is much greater than the ejecta mass (M^j) 
and the swept-up energy is still negligible. But, the solution may still be scalable from one SNR 
to another. Prom Eqs. (j29p - p3p . we also have iT=0 and iL=H=~ip/'i- This means that one SNR 
evolving within an ambient medium of density pa and temperature Ta can be scaled to another 
SNR of the density pb = /CqA"'^*^ but of the same temperature; these two SNRs have their ages 
linked by tb = ta\^^ and have the same swept-up masses and energies. 

The same scheme introduced in the previous section can also be used for embedding SNR 
seeds including the ejecta. But in this case the library of SNR templates needs to be expanded 
because the scaling is now only accurate for SNRs evolving under the same ambient temperature. 
We can tabulate a series of SNRs simulated for a temperature grid. An interpolation may be used 
to generate any needed seed for a particular ambient gas temperature. If the grid is sufficiently 
fine, then such interpolation may not even be needed. Por example, a logarithmical grid interval of 
0.02 (i.e., only 50 SNR simulations are needed to cover an order of magnitude temperature range) 
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would introduce an uncertainty of < 2% in the ejecta mass, if the template with the nearest grid 
temperature is used. Such a small variation of the ejecta mass has a negligible effect on the SNR 
inner structure. 

As indicated in §2.2, the scalability of an SNR solution also requires that its initial condition 
(i.e., SN ejecta profiles) to be scalable with respect to the surrounding medium. We find that 
such an initial condition can be set up within the uncertainty of SN ejecta models. We adopt the 
density and velocity prof i les of a post-deflagration stellar remnant of a Type la SN as proposed by 



Dwarkadas k Chevalier I (|l998l ^: 

p{r) = pse'-'/'\ v{r)=vs-, (34) 



where Ps and Vg are the corresponding values at the characteristic radius r^. The ejecta extends to 
a radius r, so that 







4:7rr^p{r)dr = Mej, (35) 
2TTr^p{r)v{rfdr = Esn- (36) 







Outside Tj is the ambient gas with an assumed uniform density pa- To make the initial condition 
scalable, we set two dimensionless parameters, 

fi = pin)/ Pa, (37) 

and 

U = ^nparf/Mej, (38) 

to be the same for all SNRs in the consideration. The four Eqs. ()35p - (|38p thus determine the four 
parameters: rs,ri,ps, and Vg. From these equations, it is also easy to show that 

2x\e'/^ -l)-2x^-x = {2,fdm)-\ (39) 



P \ 1/2 



1.5 fmhMejP J ' 

where x = rg/vi and /? = 24x^(e^/^ — 1) — (24^;^ -|- 12x + 4 + x~^). Thus x and Vg depend only 
on the assumed constants, fi and fm- Similarly, the ratio, Pg/pa = fi^^^^~^, is again the same for 
various ambient densities. Thus, we can get any desirable SNR from a pre-simulated template with 
the above scalable initial condition. 

To make the initial free expansion a good approximation to be described by Eq. (|34p . we need 
to have fm much less than one (e.g., 10~^ in our examinations; no significant difference is found if 
/m=10~^). The parameter fi (adopted to be 10) determines the shape of the initial ejecta profile; a 
larger fi (which would result in a larger r^), for example, and would give a flatter ejecta profile (i.e., 
more ejecta mass is distributed near rj). But different choices of fi (between 1 and 100) produce 
negligible effects. The same method can also be u sed to produce other form s of scalable initial 



ejecta profiles, e.g., a power law p{r) oc r " (e.g., iTruelove &: McKed Il999l ). or even a uniform 
distribution. 
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3.5. SNR Reverse Shock 

With the scalable initial condition for the SN ejecta, we can further study how the evolution of 
an SNR reverse shock depends on various physical parameters. We demonstrate this by studying 
the return time of the reverse shock (t^, i.e., when it reaches the center). Specifically, we examine 
the relation between t/j and the SN ejecta mass M^j. In general, this relation cannot be determined 
in a pure analytical form, but can be easily identified in simulations. 



Ferreira &: de Jageij (j2008l ) have examined the relation based on a series of simulations, in 
which the SN ejecta is initially distributed uniformly within a radius of 0.1 pc and has a radial 
velocity increasing l inearl y outward. They show oc M^j'^, in contrast to oc M^^^^ pred icted by 



Truelove &: McKed (|l999l ) from a simple dimensional analysis. iFerreira Sz de Jageij (|2008l ) suspect 
that this deviation may be caused by the non-zero ambient temperature assumed for the SNRs in 
their simulations. However, we find that the deviation is most likely due to their choice of the initial 
ejecta distribution, which is not scalable. Their initial condition for their — M^j examination 
requires iE=0, ip=0, and iR=0 (due to the specific choice of the initial ejecta radius), hence iM=0. 
Therefore, each of their simulations is specific to a particular choice of Mgj and is not scalable to 
different M^j value. 



Using the scalable initial condition introduced in Eqs. (j34p - (j36p . the simulations become scal- 
able. Given iE=0 and ip=0, it is easy to show it = ^im (i-e., tji oc M^/^). Based on four testing 
simulations with different M^j we identify their return times. The simulated relation of tn versus 
Mf>j is shown in Fig. [H This result is exactly the same as the expected from the scaling relation. 
The relations of tji versus Esn and pa can be obtained similarly. Finally, we have 



tR 



loM — 

rur 



-1/3 



M, 



1.4 Mm 



5/6 



■py yi.T:xvxQ/ ylO^^ergS 

p — i.ui j^yj gem ^. The same scaling relation was also obtained for the revere shock to 
reach the mantle of a core-collapse supernova exploded in a uniform medium (jRevnolds &: Chevalier 
19841 ). 



1/2 



year. 



(41) 



where m„ = 1.67 x 10' 
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4. Discussion 

We have described how a hydrodynamic solution or simulation may be scalable and how the 
scalability may be used to find out the underlying dependence on various physical parameters. In 
particular, we have demonstrated how to apply the scaling method to adaptively generate SNR 
seeds in large-scale 3D simulations of the ISM. We have also discussed how an assumed initial 
condition may affect the scalability, and specifically how the initial ejecta mass and its distribution 
are related to the return time of the SNR reverse shock. 

Potentially, the scalability can be applied to a broad range of topics. In the applications that 
we have discussed, the ambient medium is assumed to be uniform on the relevant scales. But the 



1 



Fig. 1.- 



The relation between the reverse shock return time tn and the ejecta mass Mgj for M^j 



0.7, 1.4, 2.8, and 5.6 M©. The soh d Une denotes the expected relation tn gc M^/^, while the dotted 
line the relation tR oc M~^/^ from Ferreira &: de Jager ( 2008 ). 



scaling method is still valid for any ambient medium with scalable profiles such as a power law 
density profile (e.g., p oc generated previously by a stellar wind). The medium may also be 
clumpy. As long as the inhomogeneity does not significantly affect the overall dynamics, which is 
normally true in the early stage of SNR evolution, the scaling method may still be applicable. 

We have focused on SNRs in the hot tenuous medium for simplicity and for the need of our 
practical research projects on galactic bugles. But such SNRs are not limited to those from Type 
la SNe. Most of SNRs from core-collapsed SNe may also evolve in hot gas within superbubbles, 
because massi ve stars are born m ostly in OB associations. This kind of SNRs are typically difficult 
to detect (e.g.. |Jiang et al.ll2007l ). except for situations in which radiat ion from pulsar wind nebulae 
dominates (e.g., Crab Nebula and SNR 054.1+0.3; ILu et al.l |2002| ) . Luminous SNRs that are 
dominated by shock-heated hot gas typically originate from run-away stars and happen to be in 
a relatively dense ambient medium. Such SNRs probably represent a minority of the entire SNR 
population. In the late evolution of such an SN R, the cooling becoines iinportant. Even in th is 



case, the scaling method may still be useful (§ 2.2; ISgrdll972l : IChevalierlll974l : IShelton et al. Ill999l ) 



One may also find useful applications of the scheme that we have developed to adaptively 
generate SNR seeds and to embed them into large-scale 3D simulations of the ISM. In particular, 
existing simulations of the structure formation in the universe typically use various recipes to model 
the subgrid astrophysical processes. Such recipes are often hardly calibrated with any observations 
and/or are implemented in over-simplistic ways, constrained by the limited dynamic ranges avail- 
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able in these simulations. We believe that this problem may be circumvented by the application 
of a scheme similar to ours, which allows for a more realistic modeling of the subgrid evolution of 
important processes (e.g., individual SNRs, superbubbles around massive stellar clusters, super- 
wind bubbles around galaxies, and feedback from active galactic nuclei into the intragroup/cluster 
medium). Bridging such subgrid evolution to the global hydrodynamics of the structure formation 
is badly needed to bring the simulations closer to the reality. 

We thank Bill Mathews, R. A. Chevalier, and R. Shelton for useful comments on an early 
draft of this paper. The work is supported by NASA grants NNX06AI18G and TM7-8005X (via 
SAO/CXC). 
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